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Abstract. We investigate how large lepton asymmetries affect the evolution of the 
early universe at times before big bang nucleosynthesis and in particular how they 
influence the relic density of WIMP dark matter. In comparison to the standard 
calculation of the relic WIMP abundance we find a decrease, depending on the lepton 
flavour asymmetry. We find an effect of up to 20 per cent for lepton flavour asymmetries 

h = o(o.i). 
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1. Introduction 

A wide variety of observations point to the existence of dark matter. Several different 
candidate particles have been proposed with masses varying from 10 _6 eV for Axions 
to 10 16 GeV WIMPzillas, as well as much more massive composite objects. Maybe the 
best-motivated candidate is the Weakly Interacting Massive Particle (WIMP) with a 
mass m x = 0(GeV — TeV). One reason is that the so-called hierarchy problem of the 
standard model suggests the existence of an additional particle with a mass around the 
weak scale. Such a weak scale particle, if produced thermally, would have an abundance 
today, similar to the measured dark matter density pp. The relic abundance after 
statisticaQfreeze-ont depends on the particle mass, the freeze-out temperature, and the 
annihilation cross section. The freeze-out temperature is related to the particle mass, 
Tf ~ m x /25. Below this temperature, WIMP production ceases and the WIMPS react 
only kinetically with the remaining Standard-Model particles. The kinetic freeze-out is 
delayed and happens at a temperature of a few MeV [2]. 

| We distinguish here between statistical (chemical) freeze-out of number-changing reactions, and 
kinetic (thermal) freeze-out. 
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Assuming that all dark matter is made of WIMPS, without asymmetry between 
WIMPs and anti-WIMPs, and that there are no complications such as co-annihilation, 
one can use the WMAP 7 yr constraint on the relative relic abundance [3] 



to fix the annihilation cross-section as a function of WIMP mass m x . Several 
observations of candidate dark- matter annihilation products, like e , (anti-)protons, 
(anti-)deuterons, (anti-) neutrinos and photons constrain the dark matter cross-sections. 
For a recent overview see for example [4]. The standard calculations, as described above, 
assume no asymmetry between neutrinos and anti-neutrinos. 

In this paper, we investigate how a neutrino asymmetry affects the WIMP 
abundance. Dropping the assumption of equal numbers of neutrinos and anti-neutrinos 
leads to additional contribution to the total energy density in the early universe and 
thus potentially alters the WIMP relic abundance. This slight modification of the 
Standard Model, which might emerge from a fundamental model of neutrino masses 
and interactions, is well within observational limits and might play a crucial role in 
determining the relation between the relic abundance of a WIMP and its annihilation 
cross- sect ion. 

Very little is known about neutrino (flavour) asymmetries, due to the lack of any 
direct observation of cosmological neutrinos and the large number of different theories 
of lepto- and baryogenesis. Thus, the cosmic background of neutrinos may hide a 
large neutrino-anti-neutrino asymmetry, orders of magnitude larger than the baryon 
asymmetry j5], b = nb ~ Ub ~ 0(1O~ 10 ). Here and denote the number densities of 
baryons and anti-baryons, respectively, while s is the entropy density. 

All experimental bounds on the neutrino (flavour) asymmetries are inferred 
indirectly from measurements of the primordial abundances of light elements, the cosmic 
microwave background (CMB), and large scale structure (LSS) and combinations of 
them. 

When the universe was roughly 0.1 seconds old and had a temperature of a few 
MeV, the neutrinos were in thermal equilibrium with the photons and e ± via processes 
like e + + e~ -B- V{ + with i = e,fi,r. Subsequent cooling of the universe allowed the 
neutrinos to kinetically decouple from the plasma. As muon and tau neutrinos interact 
via neutral currents only, they decouple first. The electron neutrinos interact further 
via charged-current interactions with the ambient electrons, neutrons and protons. For 
this reason, an asymmetry in the electron neutrinos influences the neutron-to-proton 
ratio at weak-interaction- freeze-out, and thereby the 4 He-abundance. Additionally any 
lepton flavour asymmetry could be detected through its additional contribution to the 
total energy density of the universe, often expressed as the number of extra effective 
relativistic degrees-of-freedom of the plasma 



tt nM h 2 = 0.1109 ±0.0056 



(1) 



f=e,n,r L 

where JJ^l/j: £1X6 the chemical potentials of the three neutrino flavour. 




(2) 
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To be more specific, the neutrino flavour asymmetry influences the 4 He abundance 
in two ways: Increasing fi Ue would lead to a smaller n/p fraction at the freeze out 
of the weak interaction rates and hence to a smaller 4 He abundance. But this effect 
can be compensated by increasing simultaneously \fi v | and/or \fi Ur \- This would raise 
the expansion rate and thus lead to a higher freeze out temperature for the weak 
interactions and increase the 4 He abundance. Note, that the chemical potentials in 
AiVeff enter quadratically and thus a measurment of AiV e g- does not constrain the sign of 
the individual flavour asymmetries. Thus even a vanishing net lepton asymmetry might 
give rise to interesting and non- vanishing lepton- flavour asymmetries, i.e. I = ^2f If = 0. 
The neutrino flavour asymmetries can be played against each other to match the 
observational bounds on light element abundances [TJ. 

After the epoch of Big Bang Nucleosynthesis (BBN), the total energy density of the 
universe consisted of photons and the three types of relativistic neutrinos. Altering that 
energy density with a large AiV eff , e.g. large fi Uf /T, would affect the CMB spectrum 
by changing the redshift of matter radiation equality and by changing the amount of 
anisotropic stress [6] . The power spectrum of density fluctuations on small scales would 
be suppressed, leading to observable effects in the large scale structure [6]. 

The BBN bounds on the pseudo-chemical potentials £f := \i Vj jT and on AiV e g 

depend on assumptions regarding the efficiency of neutrino flavour equilibration via 

neutrino oscillations. Before the onset of neutrino flavour oscillations (i.e. at T > T osc ), 

each individual lepton flavour, 

n f - n f + n vj - n Vj 
l f = , [6) 

is conserved. After the onset of neutrino flavour oscillations, only the total lepton 
number I = Ylf= e firh remains conserved, while the individual flavour asymmetries 
vary. It might be that //-oscillations ensure the full equilibration of three initially 
different flavour asymmetries [131 O C3 ESj , such that the asymmetry in the electron 
neutrinos measured today, would be the same as in the muon and tau type. In [T3l IT5] . 
it is shown that even large primordial asymmetries with fi Ue /T = 0, \i v JT = —0.1, 
and fiv r /T = 0.1 equilibrate at roughly T ~ 4 MeV. Thus, one can assume that, at 
least for [i Uf /T = 0(0.1), neutrino flavour equilibration happens at T ~ 10 MeV, well 
before the onset of BBN |16j . As a consequence LSS, CMB, and BBN are blind to 
the earlier differences in the individual lepton flavour. For this case (/ = J2fh = 3^e), 
the observation of primordial abundances and WMAP data constrain all three pseudo- 
chemical potentials: < 0.023 ± 0.041 [5], assuming the lepton asymmetry to be 
constant between nucleosynthesis and photon decoupling. The recently reported He 4 
abundance by the ACT collaboration [9] and the WMAP7 data release lead to a 
significantly larger range —0.14 < £y < 0.12 [TO] . 

Considering only partial equilibration of three initially different flavour asymmetries 
via oscillations (i.e. fi ue ^ fi Ufl ^ Av), weakens the bounds on £/. Assuming \i Ve \i Vj 
and an effective number of neutrinos N e fi = 3.3+qq while neglecting flavour equilibration, 
leads to the bound < 2.34 (here / = /x, r) and a total lepton asymmetry of |/| < 5 
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[5] . We conclude, large lepton asymmetries and large lepton flavour asymmetries before 
the onset of neutrino oscillations are compatible with current observations. 

The consequences of lepton asymmetries on the cosmic QCD transition at T qc d ~ 
200 MeV are discussed in [29]. There it is shown that large neutrino asymmetries 
influence significantly the dynamics of the QCD transition. As the electric charge, 
baryon and lepton flavour are preserved at the same time, a non-negligible lepton 
asymmetry can induce a non-vanishing baryon-chemical potential, which in turn 
influences the pseudo-critical temperature or even turns a crossover transition into a 
first order transition. 

In this work we investigate an even earlier event: statistical freeze-out of the 
abundance of weakly interacting dark matter. We show that large lepton (flavour-) 
asymmetries (|/| > 0.01) give rise to sizeable effects on the relic abundance of WIMP 
dark matter. 

In the following chapter we we will give an (incomplete) overview of arguments how 
large lepton (flavour) asymmetries can be generated. In the third chapter we show then, 
how large lepton (flavour) asymmetries influences the thermodynamic description of the 
early universe. We give analytic estimations and numerical results and show how these 
effect the relative relic abundance of WIMP dark matter in the fourth chapter before 
we conclude our work in chapter 5. Throughout the article we set c = H = ks = 1. 

2. A short review of models with large lepton (flavour) asymmetry 

Several ideas of baryo- and leptogenesis have been proposed that lead to a universe with 
large lepton flavour asymmetries but small baryon asymmetry. For example, in grand 
unified theories asymmetries |/| ^> b can be realized based on gauge invariant initial 
particle asymmetries [18] . Another mechanism is provided by supersymmetric theories, 
c.f. [25], which have the intriguing feature of flat directions. Made up of squark or slepton 
fields, these directions can carry baryon and/or lepton number. During cosmological 
inflation these squark and sleptons are free to fluctuate and to form scalar condensates, 
carrying baryon and/or lepton number, and to release these charges by decaying to 
Standard- Model particles. This is called the Affleck-Dine mechanism [191 12Q], and in 
principle can give rise to large asymmetries in baryons or leptons. 

In the early universe, large asymmetries have to be produced after an inflationary 
phase, otherwise they diluted away by the inflationary expansion and washed out 
by the huge increase in entropy density during post-inflationary reheating. Another 
possibility to dilute a large lepton (flavour) asymmetry are sphaleron transitions at 
high temperatures. These transitions violate b and I, but preserve b — I. Sphalerons 
are the reason for suggesting 0(b) ~ 0(1). Note that this does not necessarily mean 
0(b) ~ 0(lf). For instance, two flavours could be orders of magnitude larger, if their 
sum almost cancels, U ~ — lj for lepton flavours i ^ j. 

However, sphaleron transitions might never be in equilibrium and the argument in 
favor of I = 0{b) would disappear. It was shown that, if the total asymmetry Ylfh 
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is larger than a critical value l c ~ 10~ 2 then electroweak symmetry is never restored 
and sphalerons are suppressed for all times [22] 124"] . A large lepton (flavour) asymmetry 
would survive until today. It is remarkable that recent bounds on the electron-neutrino 
asymmetry from WMAP-7yr alone and combined with the ACT data are in the vicinity 
of this limit. 

In [2TJ the authors combine the Affleck-Dine mechanism with the suppressed 
sphalerons. They consider a model with a sneutrino condensate to generate large l- 
asymmetries, I > l c , but no baryon asymmetry. In their specific model they take 
the minimal supersymmetric Standard-Model with three right-handed neutrino singlet 
superfields. For a mass of the lightest neutrino of O(10~ A eV) and a preferred mass of the 
lightest right-handed neutrino of C(TeV), their model leads to asymmetries I ~ 10 -2 
to 1. In [23j [24] it is shown how the electroweak symmetry no n- restoration due to a 
large /-asymmetry works with a small baryon asymmetry. How to regulate the baryon 
asymmetry in Affleck-Dine mechanisms is discussed in [27]. An interesting model is also 
discussed in [28], where l e = — Z M and l T = 0. 

3. Thermodynamics in the early universe 

It is an excellent approximation to assume that entropy is conserved in the early 
Universe, even through several phase transitions [34]. Even if the QCD transition were 
of first order, the amount of entropy produced would be tiny [331 134"] . It is convenient 
to work with specific densities of conserved quantities. The specific baryon asymmetry 
b is related to the baryon number density n#, via b = ub/s, where s denotes the 
entropy density. Constraints from BBN and CMB are often formulated in terms of 
f]B = Ub/ti^ = bs/n T Note that ^(tbbn) 7^ ^(Wb), while b is constant. 

Before neutrino decoupling, the radiation fluid of the early Universe is a 
tightly coupled plasma in thermal and statistical equilibrium with respect to strong, 
electromagnetic and weak forces. It is characterized by a single temperature T and 
chemical potentials fii for each particle species (with mass rrii) of the Standard-Model 
{i = quarks/hadrons, leptons and gauge bosons and we count particles and antiparticles 
separately). The chemical potentials of anti-particles m = —fii, as long as particles are 
relativistic (T > mj/3). (The chemical potentials of the non-relativistic particles (e.g. 
baryons) are irrelevant, only because they are apparently too small for these species 
to contribute appreciably to the energy density, pressure and entropy of the Universe 
during the epochs in question.) 

In this work we consider temperatures between the electroweak and the QCD 
transition, T cw > T > T qcd , where T cw ~ 200 GeV and T qcd ~ 200 MeV. For WIMP 
masses between about 5 GeV and 4 TeV, it is within this temperature range that WIMP 
annihilation freezes-out and the WIMP number density is determined. 

The equilibrium distribution function fi of a particle i and its anti-particle i is given 
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for bosons (— ) and fermions (+) by the Bose-Einstein and Fermi-Dirac distributions, 

exp^jP- =f 1 exp^F 1 =p 1 

where E denotes the energy of the particles. We are interested in the difference of the 
number densities of particles and their anti-particles, i.e. their net number density, 



fi(E) = e=z ~ i » fi{E)= » +IH _^ (4) 



ni = f- 2 / £(£ 2 - m 2 ) 1/2 [/;(£) - ME)] dE, (5) 

2?r J m . 

where gi(— gi) denotes the number of helicity degrees of freedom of a particle species 
(e.g. g e — 2). For ultra-relativistic fermions the net number density simply becomes 



^ = —HiT 
o 



m 2 



(6) 



All particles and anti-particles contribute to the energy density of the Universe, 
which is given by 

^lf £2(i?2 - m 5 1/2/j(i!)dE ' (7) 



3=1,1 



where the index j can stand for any particle or anti-particle species. 

In equilibrium, electromagnetic and weak interactions provide a set of relations 
between the chemical potentials. It follows that all gauge bosons have vanishing chemical 
potentials, and for relativistic fermions 

Ik = ( 8 ) 

Hu + Hf = Hd + Hv r (9) 

A*u = At c = fit, (10) 

fid = fis = fi>b- (11) 

The indices u, d, c, s, t and b represent the up, down, charm, strange, top and bottom 
quarks. Lepton flavour (e, fi, or r) is denoted by /. 

This leaves us with five independent chemical potentials (i.e. three for the neutrino 
flavours and two for the up and down quarks) plus the temperature. These six variables 
are uniquely determined by five conservation laws and the Friedman equations. The five 
conserved quantities in the energy range of interest (between T cw and T qcd ) are electric 
charge (q), baryon number and three lepton flavour numbers. 

sq = - Ui + \ Yl Ui ~\ Yl Ui ( 12 ) 

i=e,fj,,r i=u,c,t i=d,s,b 

sb = - Yl Ui ( 13 ) 

i=u,d,c,s,t,b 

slf = n f + n Vf , f = e,fi,r (14) 

Lepton flavour is conserved, since the timescale for ^-oscillations is much larger then the 
Hubble time for T > few MeV. Additionally we assume that these global conservation 
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laws hold locally, i.e. there are no electric currents and no baryon or lepton (flavour) 
diffusion. 

We assume a charge neutral universe (q = 0) and fix the baryon asymmetry to 
b = 9 x 10" 11 [5]. While these are well established assumptions, measurements of 
the lepton flavour asymmetries in the early Universe are not available. We thus keep 
them as free parameters. In [29] it is shown in detail how to extract all chemical 
potentials fa(T; b, {If}), including details on the numerical method used to solve the 
five conservation equations function of T. 



3.1. Effective relativistic degrees of freedom 



Let us now take a closer look at the contribution of lepton flavour asymmetries on the 
effective relativistic degrees of freedom contributing to the total energy density 

30 



0*(T, {//*}) = 



"tt 2 T 4 



t{T,{fa})- 



(15) 



Together with the solution fa = fa(T; b, l e , l^, l T ) we find 

g* = g*(T;b, Z e ,Z M ,i T ). (16) 

In the ultra- relativistic case (m« = 0) and for vanishing baryon and lepton asymmetry 
(i.e. vanishing chemical potentials) we recover 



0.(T,{O})= 9i + l E 9i- (17) 

i=bosons i=fermions 

Large lepton flavour asymmetries generically lead to large chemical potentials of 
all fermion species. For ultra-relativistic fermions, (J7J) can be solved exactly. Assuming 
rrii = and fa = —fa leads for the energy density of a fermion species to 



30 v 



7 15 (JM\ 2 , 15 fJM_V 

8 4 WT/ 8 \irTJ 



This leads to an increase of g*, due to non- vanishing chemical potentials, 
9*(T, {fa}) = g m (T, 0) + Ag*(T, {fa}), 

with 



Ag.(T, {fa}) = Y,9* 



15 / fa y 
4 \ttTJ 



15 



8 VttT 



(19) 



(20) 



Any nonzero Ag* would therefore increase the total energy density and thus the Hubble 
expansion rate. This increased expansion rate would alter the relic abundances of 
light elements, and so must be checked against observations. The increased abundance 
of relativistic species is normally expressed as an increase of the effective number of 
neutrinos N VeS . iFrom the LEP measurement of the decay width of the Z° boson, we 
expect three active neutrino species with masses well below the electroweak scale, thus 



N, 



30 x ^ 



(21) 
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We also have to take a closer look at the /i- dependence of s/T 3 . In kinetic theory, 
including quantum statistics, the entropy density in terms of the distribution function 

/is m 

d 3 p 



J [/In/ T (1 ± /)ln(l ± /)] 



(22) 



where upper and lower signs refer to boson and fermion statistics, respectively. 

For ultra-relativistic particles, vanishing chemical potentials and all particle species 
at temperature T, we recover the well known result 

2vr 2 



s(T) 



m,ft=0 



15 



T 9*. 



Taking chemical potentials into account leads to extra contributions 



s(T, ^ 



m=0 



2T 
45 

2tt_ 2 
45 



T 



9* 



15 



j=fermions 



. Pi V 

9 * 



T 3 (g* + Ag s 



(23) 



(24) 



(25) 



Note that Ag St ^ Ag* and, in contrast to the energy density, terms quartic in the 
chemical potentials do not show up in the entropy density. When T ~ m,, we have to 
calculate the entropy density numerically. 



3.2. Analytic estimates and numerical results 

For an analytic estimate of the effect of lepton flavour asymmetries, we neglect all masses 
of quarks and leptons and assume that all lepton flavour asymmetries are small enough 
to justify Hi/{jiT) 1. Further we assume m w /3 > T > m b /3, thus g* = 345/42. 
From the conservation of charge, baryon number and lepton flavour we find 





23_7T 2 

23vr 2 



-T s b 



6 



T% 



^T 2 (4/i u - 3/i d 

-r 2 (2/i u + 3^ d ) 



f^e fJ'n 

+ 0(bfilfi 3 ), 



Solving this set of equations results in 



TTT 

H_ 
irT 

TTT 



7T —0 / 

2 3 

7T [26 + /] , 



1 5, 23, 

7T -b 1 H It 

6 9 3 ; 



1, 10, 23, 

7T — b-\ 1 H If 

3 9 3 * 



(26) 
(27) 
(28) 

(29) 
(30) 
(31) 
(32) 



Lepton (flavour) densities large compared to the baryon density affect not only the 
number densities of leptonic species, but also those of quarks. This can lead to increase 
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1 96^ 

Ag*(T, 6 = 0,/ / = 1/3) = -^n 2 l 2 « 6.2 x 10 3 / 2 . (34) 



of the effective degrees of freedom in the energy density and the entropy density. Note 
that all relativistic particle species equipped with baryon or lepton number contribute 
to this effect, as long as they are in statistical equilibrium. 

3.2.1. Flavour symmetric lepton asymmetry Let us first assume that all lepton flavour 
numbers are the same, l e = Z M = l T = 1/3. We also assume 6 <C \l\ <C I. Thus we can 
put 6 = 0. This results in 

^ = ~-l, ^ = -/, ^ = 2vrZ, ^ = — I, (33) 
ttT 3 ' ttT ' txT ' txT 3 y ' 

and allows us to estimate the change in the effective degrees of freedom, 
\l 

Thus Ag^/g^ « 760/ 2 , which we assumed to be small for the purpose of the analytic 
approximation, i.e. I < 10 -2 . 

For lepton asymmetries / > 10~ 2 we rely on a numerical solution of the five 
conservation equations and include all particles from the Standard-Model of particle 
physics with their measured physical masses (the unknown masses of the Higgs and the 
neutrinos are irrelevant in the regime of interest). We solved the equations f fT2|) to (fT4l) 
using (jSJ) and ( 1221) with the method described in detail in [29] . 

The numerical results for the effective helicity degrees of are shown in Fig. [TJ We 
found that an asymmetry If = 0.01 leads to a small deviation from the standard case 
with b = I = If = 0. If we apply the experimentally given upper bound for the electron 
neutrino asymmetry to all flavour, we found for If = 0.1 approximately additional 50 
degrees of freedom in the early universe between 1 < T < 50 GeV. 

3.2.2. Flavour asymmetric lepton asymmetry Let us now have a closer look at scenarios 
in which at least one of the three flavour lepton numbers satisfies \lf \ S> b, but we restrict 
to | lf \ < 1 for all flavour. For simplicity we can put 6 = 0. The first interesting situation 
is that one flavour asymmetry dominates, say If and the other flavour asymmetries 
vanish. In that case we would find that the quark chemical potentials are affected: 



ttT 


= 

3 


(35) 


irT 




(36) 


T 


57T 


(37) 


T 


64tt 
9 


(38) 


T 


10tt 

~ 9 r ' 


(39) 


/V 
T 


797T 

- 9 l - 


(40) 



The numerical results for this situation are presented in Fig. |2j We see again a tiny 
deviation from the standard case for l T = 0.01. For l T = 0.1 there would be around 10 
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T [GeV] 



100 



Figure 1. The numerical solution for the flavour symmetric case I = 3lf. The effective 
degrees of freedom of all particles in statistical equilibrium versus the temperature in 
GeV on logarithmic scale. The black line corresponds to standard case, where lepton 
asymmetries are neglected. The blue line shows the influence of If = 0.01, the green 
If = 0.05 and the red If = 0.1. 



more degrees of freedom. The effect is smaller compared to symmetric case since the 
total I is smaller. 

Another interesting scenario is a vanishing total lepton asymmetry / = 0, but 
= —l T 7^ 0. In this case quark chemical potentials would not be affected: 

fJ'Uf 237T , 



A*d — A*u — A*e — fa 



0. % 

' TFT 



«T .1 ■ — < 41) 
We find Ag* = 5(237r) 2 /^ « 2.6 x 10 4 / 2 .. As we assumed for the analytic approximation 
that the modification is small, its regime of validity is limited to \l T \ < 10~ 2 . The 
numerical results for vanishing I, but non-vanishing lepton flavour asymmetry are 
presented in Fig. [3j In the calculations for the degrees of freedom the sign of a possible 
asymmetry does not play any role, since they enter squared. For 1^ = —l T = 0.1 
additional 20 degrees of freedom appear. What makes this scenario the most interesting, 
is the possibility for even larger asymmetries. If we assume Z M = — l T = 1 we find more 
then 600 additional degrees of freedom. 



4. Effect on decoupling of WIMP dark matter 

The WIMP is a very well- motivated candidate to solve the dark matter puzzle [32]. In 
the following we assume the WIMP to be the single component dark matter particle, 
without asymmetry between WIMPs and anti-WIMP. In the hot early universe these 
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T [GeV] 



Figure 2. The numerical solution for the flavour asymmetry I = l T . The effective 
degrees of freedom of all particles in statistical equilibrium versus the temperature in 
GeV on logarithmic scale. The black line corresponds to standard case, where lepton 
flavour asymmetries are neglected. The green line shows the influence of l T = 0.01, the 
blue stands for l T = 0.05 and the red for l T = 0.1. 




1 10 100 

T [GeV] 



Figure 3. The numerical solution for the flavour asymmetry 1^ = —l r and I = l e = 0. 
The effective degrees of freedom of all particles in statistical equilibrium versus the 
temperature in GeV on logarithmic scale. The black line corresponds to standard case, 
where lepton flavour asymmetries are neglected. The green line shows the influence of 
l T = 0.01, the blue stands for l T = 0.05 and the red for l T = 0.1. 
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particles with masses typically between 10 and 1000 GeV are in thermal and statistical 
equilibrium with the radiation content. Their statistical freeze-out, when WIMPs 
decouple statistically from the radiation plasma, happens at Xf ~ m/25, leading to 
a mass dependent interval of 0.4 GeV < T fo <40 GeV (see e.g. [2]). In the following we 
want to investigate the effect of changing the Standard Model boundary condition by 
allowing large neutrino asymmetries on the WIMP freeze-out. 

To calculate the relic density, one assumes annihilations of the WIMPs XX 
with a typical weak interaction cross section o oc G\. The corresponding Boltzmann 
equation then leads to an equation for their net particle density 

n + 3Hn = - (aM) (n 2 - n 2 q ), (42) 

with the Hubble parameter H, the total annihilation cross section a and the relative 
velocity of the annihilating particles v. The index "eq" indicates the assumption of 
thermal and statistical equilibrium. Since the WIMP is non relativistic at the time 
of decoupling, the annihilation cross section can be approximated for most 
(av) = a + b/x + 0(x~ 2 ), where the numbers a and b describe s- and p-wave annihilation 
and x = m/T, for more details see e.g. (32]. It is common to introduce the specific 
WIMP abundance Y = n/s and to assume the conservation of entropy: It is convenient 
to analyze the freeze-out process in a comoving volume, thus we introduce the specific 
WIMP abundance Y = n/s and Y eq = n eq / s and the ratio x = m/T: 

ay = ^» y2 _ (43) 

ax 6c z s Hx 

In radiation domination 3c 2 = 1 is a good approximation. As freeze-out happens when 
the rate of annihilations drops below the Hubble rate, we find V^ ~ Yf ~ [H / ((crv)s)]{ . 
Thus an approximate solution for the WIMP relic abundance today can be derived as 



= Y(x -»• oo) « ( V^M I (44) 

V 9s* J f m Pl m(a/x {o + b/x] Q ) ' 

for more details see e.g. [TJ [36] . 

It is convenient to express the relative WIMP mass density f2 w i mp = p w im P /Pc with 
the critical density p c . The present relic abundance is then p w i mp = mn = msol^o, 
where the subscript zero denotes the value today: 

Q wimp h 2 = 2.7 x lO" — gL- y^ (45) 

8.5 x 10" 11 f ^\ x {o 

(46) 



GeV 2 V 9s* J f a + b/x {o 
The abundance of a WIMP particle is inverse proportionally to its annihilation cross 
section, a more strongly interacting particle stays longer in equilibrium. The dependence 
on helicity degrees of freedom is apparent in the denominator, but there is also an 
implicit dependence in Xf Q ~ c — ln(g*), where the constant c includes a logarithmic 
dependence on the wimp mass, its cross section and its helicity degrees of freedom [2]. 
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As long as Ag* -C g*, which is the case for If <C 1, we find 
AQ wimp If 1 \ Ag* Ags* 




(47) 



where £ = 0,1 when the s- or p-wave annihilation channel dominates. This can be further 
simplified using the expressions obtained for the ultrarelativistic degrees of freedom. In 
this case the contributions from /i 2 terms in the first two terms cancel out, 



However, this will not be the case close to mass thresholds. For the flavour symmetric 
case, with Xf Q ~ 25 and assuming p-wave annihilation, we find Af2 w i mp /f2 wim p ~ — 30/ 2 , 
i.e. a reduction of the WIMP relic density by 30% (7%) for I = 0.1(0.05). However, for 
/ > 10 -2 higher order contributions of the statistical potentials are important and the 
interesting regime must be analyzed numerically. 

For the numerical analysis we assume the governing effect in helicity degrees of 
freedom. The difference of the relative relic density introduced by lepton (flavour) 
asymmetries at a given freeze out temperature can be approximated 



A numerical analysis including all mass thresholds and not restricted to small chemical 
potentials is shown in figure H] for different lepton asymmetries //. We plot the ratio 
^wim P (^ J/)/^wim P (^ = = 0) as a function of WIMP freeze-out temperature. In 
the flavour symmetric case (figure H]) we find that the analytic estimate for lepton 
asymmetries overestimates the effect. We observe an effect of order 1 % for If = 0.01 
and of almost 20% for If = 0.1. 

5. Conclusion 

In this work we extend the study of cosmic lepton and lepton flavour asymmetries to 
the early Universe before the cosmic QCD transition, a subject studied previously [29] . 

The standard model of cosmology does not know any lepton or baryon number 
violating processes after the electroweak phase transition. A lepton flavour asymmetry 
produced before or during the transition would remain constant until the onset of 
neutrino flavour oscillations at T ~ O(10) MeV. Neutrino oscillations ensure that all 
lepton flavour asymmetries agree at the time of BBN, i.e. l e — 1^ = l T , and thus bounds 
from BBN and CMB observations apply to the total lepton asymmetry / only. These 
bounds are rather weak and allow |/| ^> b. 

Before the onset of neutrino oscillations, scenarios with different individual flavour 
asymmetries are even less constrained, e.g. l e = b and —l^ = l T ^> b. 

In this work we took a closer look on the the freeze-out of WIMP dark matter. This 
event is of particular interest, since observations of these relics would open a window to 




(48) 




ff*(U/^0) 9s*(l,lf = 0) 
g m {l,l f = 0) gsMJj^O) 



(49) 
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Figure 4. Comparison between the relative WIMP dark matter abundance without 
lepton flavour asymmetries QdmQ = If = 0) and several large asymmetries If =l e = 
1ft = l T , plotted versus the freeze out temperature Tf Q on a logarithmic scale. We 
observe for flavour asymmetries If = 0.05 an effect of approximately 7 percent. For 
If =0.1 this effect increases to 20 percent. 



the pre-BBN era. We have shown, that lepton asymmetries might have a sizable effect 
on the relic abundance of the WIMP dark matter. The relative relic abundance depends 
on the annihilation cross section, the WIMP mass and lepton flavour asymmetries. 

A reduction of the WIMP abundance due to large lepton (flavour) asymmetries 
happened due to an increase of the effective relativistic degrees of freedom of the 
radiation plasma. Large lepton asymmetries lead generically to large chemical potentials 
of all fermion species, which in turn contribute to the energy and entropy densities of 
the universe. We presented analytical estimates and numerical studies of this effect. 
Figures [U to [3] show how the effective relativistic degrees of freedom g* increase for 
different distributions of individual lepton flavour asymmetries. 

We demonstrated how a large If influences the relic WIMP abundance in figure 
HJ If the asymmetries are equal in in all flavour, l e — = l T and \lf\ = 0.1 we found 
a huge effect on the relative relic abundance of approximately 20 per cent. Even if 
the total lepton asymmetry is of the order of the baryon asymmetry, but individual 
flavour asymmetries are large before the onset of neutrino oscillations, a reduction of 
the WIMP relic density of up to 10% is possible. The effect presented in this work adds 
to the astrophysical and particle physics uncertainties to be accounted for in direct and 
indirect dark matter searches. Large lepton (flavour) asymmetries could be compensated 
by decreasing annihilation cross sections, in order to keep the relic abundance of a WIMP 
candidate fixed. 
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